Mon. Not. R. Astron. Soc. OOP. [TlfSOl (20111 Printed 12 January 2011 



(MN M?eX style file v2.2) 



The Massive Star Clusters in the Dwarf Merger 

ESO 185-IG13: is the Red Excess Ubiquitous in Starbursts? 



A. Adamo^*, G. Ostlin\ E. Zackrisson\ and M. Hayes^ 

^Department of Astronomy, Stockholm University, Oscar Klein Center, AlbaNova, Stockholm SE-106 91, Sweden 
^ Observatoire Astronomique de I'Universite de Geneve, 51, ch des Maillettes, CH-1290 Sauverny, Switzerland 



O 

u 

o 



> 

oo 



Accepted December 23, 2010. Received November 4, 2010 



ABSTRACT 

Wc have investigated the starburst properties of the luminous blue compact galaxy 
ESO 185-IG13. The galaxy has been imaged with the high resolution cameras onboard 
to the Hubble Space Telescope. From the UV to the IR, the data reveal a system shaped 
by hundreds of young star clusters, and fine structures, like a tidal stream and a shell. 
The presence of numerous clusters and the perturbed morphology indicate that the 
galaxy has been involved in a recent merger event. Using previous simulations of shell 
formation in galaxy mergers we constrain potential progenitors of ESO 185-IG13. The 
analysis of the star cluster population is used to investigate the properties of the 
present starburst and to date the final merger event, which has produced hundreds 
of clusters younger than 100 Myr. We have found a peak of cluster formation only 
3.5 Myr old. A large fraction of these clusters will not survive after 10-20 Myr, due 
to the "infant mortality" caused by gas expulsion. However, this sample of clusters 
represents an unique chance to investigate the youngest phases of cluster evolution. 
As already observed in the analog blue compact galaxy Haro 11, a fraction of young 
clusters are affected by a flux excess at wavelengths longer than 8000 A. Ages, masses, 
and extinctions of clusters with this NIR excess are estimated from UV and optical 
data. We discuss similarities and differences of the observed NIR excess in ESO 185- 
IG13 clusters with other cases in the literature. The cluster ages and masses are used 
to distinguish among the potential causes of the excess. We observe, as in Haro 11, 
that the use of the IR and the (commonly used) I band data results in overestimates of 
age and mass in clusters affected by the NIR excess. This has important implications 
for a number of related studies of star clusters. 

Key words: galaxies: starburst - galaxies: star clusters: individual: ESO 185-IG13, 
- galaxies: irregular 



1 INTRODUCTION 

ES0185-IG13 (hereafter, ESO 185) is a luminous {Mv = 
— 19.4) blue compact galaxy (BCG) at a distance of 80 Mpc 
(NEEq). The galaxy is now experiencing an active starburst 
episode. The relatively low exti nction and metallici ty across 
the galaxy {E{B - V) ^ 0.16, lOstlin et all (|200ll )) makes 
this irregular system very bright at blue wavebands. The 
galaxy shows morphological and dynami cal signatures of a 
recen t merger jOstlin et all 1 19991 : [Ostlinll2000l : [Ostlin et all 
I2OOII ): a tail in the north-east direction, numerous bright 
star clusters. The Fabry-Perot Ha velocity field revealed two 
distinct dynamical components, where the secondary com- 
ponent is counter-rotating and with an estimated mass an 
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order of magnitude lower than the mass of the ma in galaxy. 
The total stellar mass is of f» 7x IO^Mq (|6stlin et al. 200j). 
However the observed rotation velocity cannot support more 
than a quarter of this mass, and even after adding the con- 
tribution from the velocity dispersion the galaxy seems to 
have a lower kinematical than stellar mass, suggesting that 
it is not in dynamical equilibrium. ESO 185 has also been 
detected by IRAS at 60 /im at 0.4247±21% Jy, and has an 
integrated Hi mass of 3 x IQ^Mq. The present starburst, 
which has produced hundreds of star clusters, is probably 
the result of a minor merger between a more massive and 
evolved system and a gas rich dwarf galaxy, which has sup- 
plied the resulting system with new fuel. 

New high-resolution Hubble Space Telescope (HST) 
imaging datifl (# GO10902, PI: Ostlin) of the galaxy ab 
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lowed us to investigate the formation and evolution of the 
system, through the analysis of the cluster population. A 
composite image of the galaxy is showed in Figure[T] The top 
panel shows the galaxy in the deep WFPC2 F606W frame 
with a 3-colors inset of the central starburst region. Three 
different regions are outlined and showed in the other two 
panels and in Figure [2] The Figure [2] shows a 3-colors com- 
position of the galactic main body. A barred-like structure 
crosses the center of the galaxy and is formed by hundreds of 
star clusters and several starburst knots (cluster complexes). 
Arm paths are traced by star clusters and surrounded by 
Ha emissio n (green halo in the image). The Ha emission is 
prominent. lOstUn et al.1 (|l999l ) estimated a Ha luminosity 
of Lhci 80 X 10'^'^ W and an ionised gas mass of few times 
10* Mq. In the bottom left panel (Figure [TJ in false colors, 
we show the tidal stream or plume (both terms will be used 
in the text) extracted from the F606W image. The plume is 
approximatively 32.0" long. The contours draw the shape of 
the plume at 2 — 3o- levels from the background. It is a low 
density structure. We detect clusters (purple circles) at the 
end of filament-like structures. We have, for the first time 
in this galaxy, detected a shell structure 8.0" away from the 
starburst region (13.5" radius from the galactic center) in 
the southwest region of the galaxy. Numerous star clusters 
are observed in the region between the shell and the main 
galactic body and on the shell. 

Shells have been observed in many early-type post- 
star buJst_gaJMiBS_aa^_^^ typic a l signatures of recen t merg - 
ers (|Hernquist fc Spergell Il992l : iKoiima fc Noguchil 
Using self-gravitating syst e ms o f stars and dissipative gas 
clouds, iKoiima fc Noguchil (|1997D simulated a minor merger 
event between an elliptical system and a gas-rich satel- 
lite. They found that the gas of the "sinking satellite" was 
not only confined to the inner regions of the new formed 
system, but also distributed in the fine structures formed 
during the merging, like shells and ripples. Interestingly, 
they observed short starburst episodes, abruptly reduced 
due to the dissipative nature o f the gas. Millimeter obser- 
vations of NGC 5128 (CenA, ICharniandaris et al.l l2000f ). 
revealed the presence of molecular gas across two of the 
many shells surrounding the galaxy. Similar features were 
also observed in the Medusa merging system (NGC4194, 
lAalto et al.ll2001^. These observations confirmed the find- 
ings by Koiima fc NoguchU (|l997l ) and suggested the possi- 
bility of cluster formation in shells. Recently, from the anal- 
ysis of the globular cl usters (GCs) in the post -starburst shell 
galaxy AM 0139-655, iMavbhate et al.l (|2007l ') found some of 
the youngest GCs (~ 0.4 Gyr) located in the shell. They 
discussed the possibility that the clusters and the shell may 
have formed in the same dissipative merger event, and that 
the star formation in the galaxy was shut down after this 
last interaction. 

The formation of star clusters is usuall y enhanced in in- 
terac ting/m erging systems (the Antennae lWhitmore et al 



19991 : M51, iBastian et al] l2005l : the Bird. IVaisanen et al 



20081 amon g many others). They have also been foun d in 
tidal tails (jKnierman et ahlbOOsI : iTrancho et al.ll2007l ') and 



coUisional rings (jPellerin et aLllioiOl ) . Starburst dwarf galax- 
ies have also pr oved to be very effi cient in prod ucing massive 
youn g clusters (iBillett et al. 1|2002'; NGC 1569. lAnders" et al.l 
|2004 NGC 1709. lAnnibah et al.,,2009i : to cite a few exam- 
ples). The BCGs belong to this wide group of irregular and 
low mass galaxies. Deep imaging stu dies of the cluster p opu- 
lations in BCGs, like ESQ 338 -IG04 (|6stUn et al.ll2003h and 
Haro 1 1 (|Adamo et aLllioiOal ) , have revealed very high clus- 
ter formation rates (CFRs) during the last 30-40 Myr, and a 
peak of cluster production less than 10 Myr ago. The anal- 
ysis of the star cluster population in BCGs can be used to: 
1) reconstruct the formation history and the environmental 
properties of the host galaxy; 2) investigate the complex and 
short embedded phase, during which more than 90 % of the 
clusters will be des troyed (the so-called "infant mortality", 
iLada fc Ladall2003l ). 



The embedded phase is a transient period in the clus- 
ter evolution. The clusters form from the collapse and frag- 
mentation of giant molecular clouds (GMCs). The strong 
UV radiation, produced by the young massive stars, ion- 
ize the intracluster gas and create Hii regions around the 
clusters. However a large fraction of stars is still accret- 
ing material from their dusty disks (young stellar objects, 
YSOs) or contr acting in the pre-main sequence phase (see 
iBik et al] 120101 for an example of different components of 
a massive star- forming region in the Milky Way). During 
this partial ly embedded phase, which could last a few Myr 
(|Ladall20ld ). the gas expulsion is the main mechanism caus- 
ing the early dissolution of the clusters. NIR studies of 
very young unresolved clusters in nearby galaxies have re- 
vealed fiuxes that are impossible to reconcile with evolu- 
tionary models of stellar continuum and photoionized gas. 
In Haro 11, we uncovered a very pr ominent population o f 
clusters affected by this N IR excess (lAdamo et al.ll2010al) . 
In other starburst g alaxies (iFernandez-Ontiveros et al. 20091 : 
ICresci et all l2005l : ICabanac et all 120051 ) . the NIR excess 
found in very young clusters has been explained as due to 
a significant flux contribution from young stellar objects 
(YSOs) or from other sources of hot dust. The use of IR 
data clearly reveals a much more complex scenario than the 
one predicted by single stellar population models and more 
investigations are required to increase our understanding of 
the star formation process in clusters. 



Space Telescope, obtained at the Space Telescope Science Insti- 
tute, which is operated by the Association of Universities for Re- 
search in Astronomy, Inc., under NASA contract NAS 5-26555. 



The paper is organized as follow: In Section [5] we de- 
scribe the data reduction and the photometric analysis, in- 
cluding completeness tests. Section |3] contains a detailed de- 
scription of the spectral energy distribution (SED) analysis 
of the star clusters, and how we constrain the properties of 
clusters affected by the red excess. The final age, mass, and 
extinction distributions are shown in this section, along with 
Monte Carlo simulation tests. We discuss the possible mech- 
anisms causing the NIR excess in young clusters in Section 
13] In Section O we discuss the environmental properties of 
the host starburst galaxy as revealed by the star clusters 
population. In the last paragraph of this section, we inves- 
tigate the formation and possible evolution of the galactic 
system. The conclusions of this work are summarized in the 
last section. 
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Figure 1. Composite of tlie starburst galaxy ESO 185-IG13. On the top left panel, a saturated WFPC2 F606W image, with a 3-color 
inset of the central starburst, clearly shows the position of the shell and of the tidal stream. The galaxy has been sampled with the high 
resolution Planetary Camera, while a part of the extended plume, in the Northeast side, falls onto the WF2 chip. Three different regions 
of the galaxy are outlined and showed in the corresponding panels. Tidal stream: Faint surface brightness contours reveal the shape of 
the plume as shown in Region (a). The purple circles indicate the position of the star clusters that we detected only in the two F606W 
and F814W deepest exposures. Starburst, Region (b); it is shown in Figure[2] Shell: In Region (c), we show a particular of the tidal-shcU 
feature visible in the Southwest side of the galaxy. Several clusters are located between the central region and the shell. A few are also 
visible across the shell. See the main text for more details. 
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Figure 2. Three-color image of the central starburst: WFPC2 F814W filter in red, F606W in green and F336W in blue. A barred-like 
structure with hundreds of SCs characterizes the nuclear region. We clearly see a hint of spiral arm features. Extended Ho emission 
(transmitted by the F606W filter) surrounds the young starburst as shown by the prominent green halo. See the main text for more 
details. 



2 DATA SAMPLE AND MULTI-BAND 
PHOTOMETRY 

2.1 Observations &: Reduction 

HST multiband high-resolution imaging of ESQ 185 were 
carried out in 2007 (associated with programs # GO 10902, 
PI: G. Ostlin.). The galaxy was sampled from the FUV 
(ACS/FUOLP) to the NIR (NIC3/F160W) in order to an- 
alyze the stellar and gas components. 

The angular size of the starburst region is quite small 
(major axis about 26.0") and easily fitted into the field of 
view of the Planetary Camera (PC), the Solar Blind Channel 
(SBC), and the NIC3 camera. In a previous short exposure of 
the galaxy (GO 5479: PI M. Malkan) the plume fell outside 
the detector. In the images we got in 2007 the position of the 
galaxy was set so that the plume could be sampled by the 
WF2 camera. The same position was used for all the data 
taken with the WFPC2 instrument. We find that the plume 
is around 32.0" long, too extended to fit into the SBC field 
of view. The plume was only partially imaged in the NIC3 
frames. 

The frames were reduced, drizzled and alig n ed us - 
ing the MULTID RI ZZLE task jKoekemoer et all (|2002h : 
iFruchter fc Hookl (|2002l )') in PyRAF/STSDA^. The final 

^ STSDAS and PyRAF are products of the Space Telescope Sci- 
ence Institute, which is operated by AURA for NASA 



WFPC2 and ACS frames were rescaled to 0.025" /pixel, 
before combining. To improve the large plate scale of the 
NIC3 camera {0.2" /pixel) the galaxy was imaged with a 
double NIC-SPIRAL-DITH path in 3 dithered pointing. 
This allowed us to reach a diffraction limited sampling, 
0.067" /pixel, of the final science frame. A more detailed 
description o f the N ICMOS data reduction is given in 
lAd amo et al.1 (|2010al ). A list of the filters, total exposure 
time, zero points (ZPs), and other observational properties 
related to the reduction of the data are listed in the Table 

m 



2.2 Cluster catalogue and photometry 

The extracti on of the position of the clusters was done with 
SExtractor (|Bertin fc Arnout3ll996t ). We cross-checked the 
detections in the two deepest exposures ( WFPC2 F606W 
and F814W). Crowded regions inside, and isolated ones sur- 
rounding the st arburst, were treate d with different set of 
parameters (see lAdamo et aDl2010ah . SExtractor was also 
run on the WF2 frame to extract potential cluster candi- 
dates in the plume. We established by eye the center of the 
starburst region (a = 19 : 45 : 05.2; 5 = -54 : 16 : 01.1) 
and considered candidate detections inside a circle of ra- 
dius 15.12". For the plume we delimited the studied region 
in the WF2 frame to a box of 18.85" x 20.8". These chosen 
projected dimensions were based on the drop of the S/N 
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Table 1. i?5T observations carried out within the framework of program # GO 10902 (PI G. Ostlin). For each filter we list in the table: 
the total exposure time; the corresponding ZPs (Vega magnitude system); the aperture corrections, the number of objects detected 
with a S/N ^ 5; and the corresponding magnitude limit, a) Both the HST filter names and the abbreviated nomenclature indicated in 
bracket are used thereafter in the text, b) The frames in the WF2 chip are used for the analysis of the plume (Region (a) in Figurc[TJ. 



Instrument 


Filter" 


Camera 


Exposure time 


ZP (mag) 


a,c{mag) 


N(o- ^ 0.2) 


mag limit 


ACS 


F140LP (FUV) 


SBC 


2584 s 


20.92 


-0.54±0.05 


210 


24.8 


WFPC2 


F336W (U) 


PC 


1200 s 


19.43 


-0.42±0.05 


96 


23.4 


WFPC2 


F439W (B) 


PC 


800 s 


20.88 


-0.36±0.05 


121 


24.7 


WFPC2 


F606W (R) 


PC 


4000 s 


22.89 


-0.61±0.09 


420 


27.5 






WF2'' 


4000 s 


22.92 


-0.31±0.05 


11 


27.9 


WFPC2 


F814W (I) 


PC 


4500 s 


21.64 


-0.73±0.04 


420 


26.7 






WF2'' 


4500 s 


21.66 


-0.35±0.05 


11 


28.1 


NICMOS 


F160W (H) 


NIC3 


4992s 


21.88 


-2.45±0.37 


126 


25.9 



outside the delimited regions. A fraction of the plume falls 
on the PC camera frame, however this region is still part 
of the outskirts of the starburst area. Therefore, we identify 
plume clusters as the ones detected in the WF2 frame (see 
Figure [1]). This cross-checked catalogue contained 1265 clus- 
ter candidates of which 14 are in the plume (shown in panel 
(a) of Figure [l|. Using this catalogue we performed aper- 
ture photometry (IRAF/PHOT task) in all the available sci- 
ence frames. The aperture radius for doing photometry was 
fixed to 0.1" (similar to the full width half maximum of the 
poi nt spread function, P SF, in the UV and optical frames, 
see lAd amo et al]|2010al for details) and a local sky annulus 
of 0.05" wide at the distance of 0.125" from the center was 
used. This set of parameters was the same for all the frames 
from the UV to the IR. We estimated then a mean aperture 
correction for each frame to take into account the flux lost 
by the fixed radius used and the different PSFs. More de- 
tails regarding the ap erture correction estimates are given in 
lAdamo et"al] (|2010al ). since Haro 11 and ESO 185 were ob- 
served with ACS/SBC, NICMOS, and WFPC2 setups. The 
aperture corrections and their associated uncertainties are 
listed in the Table [T] The cluster candidates in the plume 
were only detected in the two deepest F606W and F814W 
exposures. The two F336W and F439W frames are too shal- 
low and heavily affected by charge transfer efficiency (CTE) 
problems, therefore the plume and clusters were undetected. 
We performed photometry in the two F606W and F814W 
WF2 frames, using an aperture radius of 0.2" (comparable 
with the full width half maximum of the PSF in these two 
frames) and a local sky annulus of 0.05" wide at the distance 
from the center of 0.225". 

For the WFPC2 images, we applied CTE corrections to 
the observed fluxes using a modified version of the algorithm 
by iDolphinI (|2009l ). As already described in lAdamo et al.l 
(|2010al ). we chose not to use the median background value, 
as prescribed by Dolphin's formula, but a "local" back- 
ground for each source. 



2.2.1 Final catalogue 

The final catalogue was created applying the following cri- 
teria: 

i) the output photometric error associated with each mea- 
surement was used to reject estimates with magnitude errors 
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Figure 3. Color magnitude diagrams of the cluster population. 
Top panel: the small green (grey, in a B/W printed version) as- 
terisks show the 1265 cluster candidates detected in the starburst 
region and the plume. The black filled dots represent the 420 
clusters with photometric errors in I and R, Om ^ 0.2. The pur- 
ple (grey) filled dots are the 11 clusters detected in the plume. 
Bottom: the green asterisks show the candidates which have been 
detected in R, I, and B band. The filled dots are the clusters with 
o"m ^ 0.2 in R, I, and B band. The extinction vector indicates in 
which direction the clusters move if corrected for extinction. The 
evolutionary track from ZOl are for Z=0.008 and several values 
of the mass. See main text for details. 
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(Tm. > 0.2. In Table [T] we show the corresponding magnitude 
limits and the number of objects retained in each frame after 
this selection; 

ii) the catalogue includes only objects detected in at least 
3 filters; 

iii) the final error associated with the integrated flux in 
each filter was the combined quadrature of the photometric 
error and the uncertainty on the aperture correction deter- 
mination (see Table [l|; 

iv) we corrected the fiuxes for the foreground galactic ex- 
tinction (|Schlegel et al.ll 19981 '). 

In total the final catalog contains 290 object^ with detection 
in at least 3 filters. All the detections in the plume, except 
one (which is also detected in the H band) are excluded from 
this final catalogue. Looking at the Table [1] it is possible 
to see that in the R and I bands the number of clusters 
with Gm ^ 0.2 is significantly higher (420 in the starburst 
and 11 in the plume). We show this numerous population 
in the color-magnitude diagram (CMD) on the top panel in 
Figure[31 We also include a second CMD, replacing the deep 
I band exposure with the shallower B filter, which clearly 
demonstrates the difference in the quality of the available 
data. Sources with a S/N larger of 5 in B band correspond 
to clusters more massive than lO^M© or brighter than —10.0 
mag in absolute magnitude. 
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Figure 4. Completeness fraction as function of the object mag- 
nitudes. The outputs of the F606W frame are plotted in green 
(grey), while in black we show the fraction for the F814W frame. 
See the main text for the definition of "inside" and "outside" 
regions. The red dotted straight lines show the 90 % and 50 % 
completeness limits. 



3 ANALYSIS OF THE CLUSTER PROPERTIES 



2.2.2 Completeness limits 

We performed completeness tests on the two R and I 
WFPC2 fr ames to trace the sens itivity limits reached in the 
analysis. In lAdamo et al.l (|2010al ') , we give a detailed descrip- 
tion of the method used to construct the completeness test. 
Since the R and I imaging have the same properties as the 
ones used for the Har o 11 analysis, we use d the two PSF 
samples constructed in lAdamo et al.l ll2010ah to make grids 
of fake sources with mksynth ( LarsenI 19991 *). We created 
three different grids of positions and for each of them the 
magnitudes were assigned randomly, to avoid systematics in 
the test. As function of the magnitude bin, we estimated 
the average recovered fraction of objects (allowing a shift 
smaller than 3 pixels). The number of recovered objects de- 
pended on the crowding of the region. Similarly to what we 
did for Haro 11, we defined a radius of 6.5" to distinguish be- 
tween the crowded central region and the surroundings with 
low background. In Figure [l] we show, for each filter, the 
total fraction of recovered detections, and how this fraction 
changes as function of the environment where the sources 
have been detected. The 90 % completeness is reached at 27 
mag for both filters (continuos black and green lines). Inside 
the crowded region 10 % of objects with 25 — 26 mag are al- 
ready missed and the fraction increases for fainter ones. On 
the other hand, objects located outside the starburst area 
have been easier to detect. This result confirms that low 
mass young objects and more evolved faint clusters have 
higher chances of being detected around the intense star- 
burst region, while in crowded regions, we are not able to 
efficiently de-blend low brightness sources. 



* The complete photometric catalogue is available on request 
from the authors. 



3.1 Spectral Evolutionary Synthesis Models and 
SED fit 

To estimate the a ges, masses, and e xtinct ions of the clus- 
ters we used the Zackrisson et al. I (|200ll . hereafter ZOl) 



spectral synthesis model. This model predicts the com- 
bined SED of both stars and photoionized gas. The gas 
around the star cluster, photoionized by the young massive 
stars, contributes both to the continuum and with emis- 
sion lines to t he integrated broadband fiuxes of the young 



systems (e.g. Kriiger et al.l Il995l; IZackrisson et al] 
Anders fc Fritze-AlvenslebenI l2003l: IZackrisson et al. " 



Reines et al.ll2010l : lAdamo et al.ll2010al lbl). In lAdamo et all 
( 2010bD we show that nebular emission non-negligibly affects 
the SEDs of the clusters during the first 10 — 15 Myr of clus- 
ter evolution. The contribution becomes smaller after 6 Myr 
on the blue side of the cluster spectrum, but lasts longer in 
the NIR wavebands. This is a very important property to 
keep in mind when both optical and IR fluxes are used to 
constrain the cluster properties. 

The use of the ZOl models and the - algorithm are 
thoroughly presented in lAdamo et al] (|2010al '). For the anal- 
ysis of the ESO 185 cluster population we used the sam e 
parameter set as for Haro 11: Kroupa's IMF (|Kroupall200"ll '). 
standard H ll region values for the gas , but a slightly higher 
metallicity for the gas and the stars (Z = .008). The spec- 
tral synthesis model of iMarigo et all (|2008l . hereafter M08) 
has been used to test the robustness of some of our results. 

We performed X^"fit of SEDs with detections in at least 
3 filters. The extinction was treated as a free parameter . 
We used the Calzetti attenuation law (|Calzetti et al.|[2000l ). 
which includes a more general tr eatment of all the com plex 
environments in the galaxy (see lAdamo et alj l2010al for a 
comparison between the Calzetti and the LMC extinction 
law). The cluster mass was estimated from the normalization 
factor between the best fitting model and the observed SED. 
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Figure 5. SED analysis of 4 clusters. The filled black points indicate the observed fluxes of each cluster. The integrated model fluxes are 
labelled with different symbols for each set of fit as indicated in the plot insets. The crosses in the UV-UBR fit and UV-UBRI fit show 
which detections have not been included in the fits. We plotted the spectrum of the flnal best fitting model with the colour corresponding 
to the best fit run. The integrated model fluxes (symbols) sit above the model continuum because of the contribution by emission lines. 



3.2 Constraining the SED of the clusters; with or 
without red excess 

In lAdamo et al.l (|2010al ) we discovered that several young 
clusters have a flux excess in the I, H, and Ks bands. The 
inclusion of those filters in the fit produced high residu- 
als, and a clear overestimation of the cluster age. Due to the 
similarity between Haro 11 and ESO 185 in terms of opti- 
cal luminosity, metallicity content and patchy extinction, we 
expected to find young clusters possibly showing the same 
excess. 

In the first run of x^-fits we included detections in all 
the available filters (cr ^ 0.2). Thereafter, we refer to this 
run as the UV-UBRIH fit. However, not all the objects in 
the catalogue were detected in the 6 available frames. Of 
the 290 candidates, 50 % axe detected in 3 filters, 21 % in 
4 filters, and the remaining 29 % in 5 or 6 filters. In the 
second run of x^-fits, the UV-UBRI fit, we excluded the H 
band detections. In this case, the number of fitted SEDs 
diminished from 290 to 242. Finally the third set of x^- 
fits, the UV- UBR fit, only included the blue optical and UV 
filters. In this final run, only 117 clusters could be fitted. 
The small number of clusters is due to the shallow B and U 
frames. 



Using plots as the ones shown in Figure [S] and com- 
paring the reduced Xr and corresponding Q values in the 3 
sets of fits we could distinguish whether a cluster is affected 
by some kind of excess. In general, the exclusion of a data 
point generates a smaller Xr, however the best fitting model 
should be the same if the observed fluxes are self-consistent. 
In our exercise we found that several objects didn't behave 
as expected. As already observed in many clusters in Haro 
11, there are systems which show a better fit to the UV and 
blue optical data points when the IR and I bands fluxes are 
excluded. Systematically, the excluded observed fluxes sit 
above the best fltting models. It is clear that self-consistent 
stellar and nebular evolutionary models are not enough to 
interpret the NIR luminosity properties of these clusters and 
other mechanisms nee d to be addressed (see SectionUof this 
work and Section 5 in lAd amo et al.ll2010al ). 

In Figure [5] we show a representative sample of cluster 
SEDs in ESO 185. The x^ outputs of the 3 runs are collected 
in Table [21 More objects are included in the Appendix (Fig- 
ure [T] and [2] and Table[T]). In Table[2]of the Appendix we list 
the magnitudes, in the Vega system, of all the shown cluster 
SEDs. 

Cluster #286 is an example of an "expected" regular 
cluster SED. The top left panel of Figure [5] shows that the 
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Table 2. Final outputs given by the three sets of SED fits as 
shown in Figure [5] In bold we show the final age, mass and ex- 
tinction assigned to the clusters. 



IDs 




Q 


Myr 


IO^A/q 


E(B-V) 


UV-UBRIH?it 


286 


2.5 


0.109 


9.5 


2.6 


0.0 


212 


8.4 


5.C-05 


19.5 


4.5 


0.03 


344 


5.0 


0.018 


19.5 


2.3 


0.03 


2f5 


2.0 


0.205 


3.5 


0.4 


0.19 


UV-UBRim. 


286 


3.8 


0.056 


9.5 


2.7 


0.0 


212 


12.0 


2.e-05 


19.5 


4.5 


0.03 


344 


5.0 


0.018 


19.5 


2.3 


0.03 


215 


0.9 


0.632 


3.5 


0.4 


0.19 


UV-UBR fit 


286 


7.4 


0.024 


9.5 


2.6 


0.0 


212 


1.6 


0.444 


3.5 


0.5 


0.05 


344 


1.8 


0.401 


2.5 


0.3 


0.05 


215 


1.7 


0.430 


3.5 


0.4 


0.18 



3 sets of fits were all reasonable successful. Moreover, from 
Table [2] we see that the best age, mass, and extinction do 
not change from one fit to another. On the other hand, the 
UV-UBRIH fit fails to reproduce the SED of the clusters # 
212, # 344, and # 215. Cluster # 212 has a clear excess in 
H band, and the best model produced by the UV-UBRIH fit 
cannot properly fit the U and B data points. The exclusion 
of the H band data, in the UV- UBRI fit, does not provide 
an acceptable fit (outputs in Table [21). Only the exclusion 
of the I band allows us to find a best fitting models able to 
fit the U and B data. The outputs of the UV- UBR fit, in 
Table (5] show that the cluster is younger and less massive 
than predicted by the other two sets of fits. Cluster # 344 is 
similar in properties to # 212, except that this cluster is not 
detected in the H band. Finally, cluster # 215 has a clear 
excess in the H band but not in the I band. 

As already done in the analysis of the Haro 11 clusters, 
we estimate the properties of the clusters affected by NIR 
excess using only the UV and blue optical SEDs, or in other 
words with the UV-UBRI and UV-UBR fits. From the initial 
290 analyzed objects, 197 clusters do not have any NIR fiux 
excess and have been fitted by a UV- UBRIH fit. A total of 37 
clusters have an I and an H band (if detected) excess, similar 
to clusters # 212 and # 344; and 30 objects show properties 
similar to # 215. The remaining 26 objects, mainly detected 
in only 3 filters have been excluded because of inconsistent 
fits. Among them, 4 objects are not detected in the UV-U 
bands but show an increasing brightness from the R to the 
H filters. We will discuss the fit to these four sources in the 
following section. 

In Figure [6] we show an analysis of the residuals pro- 
duced by clusters with a flux excess. We define the residual 
as Am = mmod — n^-obs, i-e. Am gives a quantitative mea- 
surement of the offset between the fiux provided by the best 
fitting model and the observed data point. When clusters 
affected by a NIR excess (I and/or H band) are fitted by 
a UV- UBRIH fit the output residuals (see the distributions 
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Figure 6. Residuals in all the filters of clusters with H band 
excess. The showed sample includes also sources with excess in 
I and H. The hatched histograms show the residuals. Am = 
"*mod — ™obsi produced by the best fitting models determined 
with the UV-UBRIH fit. The red (grey) line histograms show the 
residuals when the UV-UBRI and UV-UBR fits (depending if the 
excess is only in H or also in I band) are used. The vertical thick 
dotted lines indicates the perfect match between model and ob- 
servation. The number of objects diff'ers between plots because of 
the different detection thresholds in each band. See text for more 
details. 



shown in hatched histograms) have large scatters around 
zero. In particular, the residuals in the H band are systemat- 
ically offset versus positive values. The red histograms shows 
the re-estimated residuals produced by the UV-UBRI and 
UV- UBR fits (i.e. excluding the data points affected by red 
excess). We clearly notice that the distributions in the H and 
I bands become even wider. On the other hand, in all the 
remaining filters, we obtain smaller residuals centred around 
zero. 



3.3 Testing uncertainties on the determination of 
the cluster properties 

We described in lAd amo et al.1 (|2010al ) the method used to 
create cluster populations using Monte Carlo simulations. 
Here, we performed the same exercise using the set of fil- 
ters available for ESO 185. Combining three values of masses 
[10'',10'\10'*] Mq and E{B - V) = [0.05,0.2,0.3], we cre- 
ated for each of the 51 time steps of the model nine differ- 
ent clusters, giving a total input mock cluster population 
of 9 X 51 = 459 objects. For each of these clusters we sim- 
ulated 1001 realisations by stochastically adding offsets in 
accordance with the observed error distribution, for which 
we fitted SEDs to determine ages, masses, and extinctions. 
We compared the output age, mass and extinction distribu- 
tions with the input parameters to test the uncertainties we 
introduced in our analysis. 

In Figure [7] we briefiy summarize the result of our test. 
Of the 459 clusters, we included in the analysis only clusters 
with apparent magnitude, m_R ^ 27.5 mag. This cut is sim- 
ilar to the one performed on the real catalogue so we could 
recreate a sample with properties comparable to our data. 
We performed four different sets of fits: a fit of the whole 
SEDs, using the UV-UBRIH fit; in the second run the UV- 
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Figure 7. Results from our Monte Carlo simulations of the uncertainties in the SED modelling. We compare the input distribution to 
the output (represented by the median of 1001 realisations). Top panels: Recovered ages (black hatched histograms). The original input 
age distribution is shown with a red (grey) line. Central panels; Recovered masses are represented with the hatched histograms. Different 
colors (and line inclinations, for a B/W printed version) are given to each input mass (dotted histograms): black shows 10^, purple 10''', 
and blue 10® clusters. The recovered masses are shown with the corresponding colors. Bottom panels: Recovered extinctions are again 
plotted with hatched histograms, while the input values were drawn in dotted lines. Black (inclined lines) correspond to the initial value 
of E{B — V) = 0.05 mag, red (vertical lines) is for 0.2 mag; and green (horizontal) for 0.3 mag. The combination of filters used in the 
SED fit are indicated in the plots. 
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Figure 8. Recovered ages as function of tiic input ones, for a RIH fit. Tliis filter combination produces very poor fits for clusters witli 
age between fO Myr and f Gyr. Clusters fainter than 27.5 mag have been removed from the analysis, therefore clusters with different 
physical parameters cover varying age ranges. Each of the panels is for a different value of mass and extinction. The dots are the median 
values of the fOOl realizations, and the error bars show the quartiles of the distributions. A dashed line in each panel draws where, in 
the plane, the input and output ages are equal. 



UBRI fit where we excluded H band; two more runs with 
a UV-RI and RIH fit, since many of the final clusters were 
fitted with these 2 combinations of 3 filters. 

In the top panels we compare the input age distribution 
(red line) with the recovered ages (black hatched histogram) 
of the simulated clusters using the 4 sets of fits. 

Fitting all the available filters or excluding the H band 
(see top left panels) give quite trustable results except in the 
region of 10-100 Myr. Clusters with input ages of 8.5 Myr 
(log(age)= 6.9) are redistributed partially at slightly older 
ages, 10-12 Myr. More important, in our analysis we tend to 
underestimate the ages of clusters between 35 and 100 Myr. 
We find, in fact, that the age bins at log(age)= 7.5 to 8.0 
are emptied while two overpopulated regions are created be- 
tween log(age)= 7.1 and 7.3. This effect has previously been 
noted in other similar multi-band analyses (|Mafz Apellanij 
l2009f l. The combination of available filters for ESO 185 was 
not able to fully break the resulting age degeneracies. In 
our analysis of Haro 11, the addition of more filters, like the 
F220W (NUV fiher), F550M (narrow V-continuum) , and Ks 
(IR), proved to be less su sceptible to this effect (see Figure 
17 in lAdamo et ahlbOlOal ). 

When only UV, R, and I bands are used, we observe 
that the gap at 30-100 Myr becomes more extended cre- 
ating a single peak at 12 Myr. This happens because the 
model track from ~10 Myr to 14 Gyr moves parallel to the 
extinction vector, creating a degeneracy between age and 
extinction. This effect can be also seen on the right panel at 
the bottom, where the recovered extinctions for the UV-RI 
fit tend to be overestimated. The same mechanism moves 



input clusters with ages around 14 Gyr down to 1 Gyr. The 
use of this combination of 3 filters could introduce impor- 
tant uncertainties in the final age distribution. We observe 
in Figure [9] that the estimated cluster ages using these 3 fil- 
ters are between 3 and 25 Myr. There is, therefore, a high 
chance that those 12 Myr old clusters are in reality much 
older and less extinguished. 

Finally, for the RIH fits, we observe significant displace- 
ments of the output age distribution (the most right top 
panel), as well as in the recovered masses and extinctions. In 
Figure |5] we show the recovered ages for the RIH fit as func- 
tion of the input ones, for different combinations of masses 
and extinctions. For all combinations of masses and extinc- 
tions, we observe that input ages 6 Myr are well repro- 
duced. Ages older than 1 Gyr are well recovered in more 
than the 50 % of the cases. Why does the RIH fit fail in re- 
covering ages between 6 Myr and 1 Gyr? The evolutionary 
tracks between 6 and 1 Gyr move parallel to the extinction 
vector and make some complex loops (Figure [T^ . This be- 
haviour lies behind the displacements we see for the lowest 
extinctions, at all simulated masses (see top panels in Figure 
At the same time, in Figure [7] bottom right panel, ones 
sees that the values of the recovered extinctions, for input 
E{B — V) = 0.05 are widely spread up to output values of 
0.4 mag. These clusters contribute to the spurious peak at 
roughly 6 Myr (log(age)= 6.8). For higher input extinction 
values the situation is less dramatic. The secondary peak at 
20 Myr (log(age)= 7.3) is created both by objects of input 
age 10 Myr for which the age is overestimated and objects 
in the usual 30-100 Myr input age range. 
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Due to the large uncertainties associated with the RIH 
fit, we could have excluded objects detected in only these 
three filters from the final analysis. However, as discussed 
in the next section, the output age distribution for observed 
clusters produced by the RIH fit (see Figure (9]) has two 
peaks, one at very young ages, and a second one at very 
old ones. For the young age bin the recovered properties are 
consistent with the input ones. Clusters as old or older than 
1 Gyr have more than a 50 % chance to be consistent with 
the real ages. Therefore, we consider that the age for most 
of the clusters observed in the RIH bands can be trusted. 

In Figure [71 the central panels shows that there is a ten- 
dency in our analysis to underestimate the masses up to a 
factor of 2. This behaviour is related to the recovered extinc- 
tions which are also underestimated (bottom panels). That 
the mass offset increases with input mass is due to the fact 
that we only include clusters with a simulated magnitude 
brighter than ma = 27.5, and to enter the sample at high 
ages they need to be massive. Clearly, the mass and extinc- 
tion uncertainties increase when only 3 filters are used. 



3.4 Results: From the embedded phases to the 
old cluster population 

The final distributions of age, mass, and extinction for the 
264 clusters are shown in Figure (9] In total, the properties of 
~50 % of the sources have been estimated from detections 
in only 3 filters: 43 clusters have detections in the R, I, 
and H bands only; 91 ones are observed in FUV, R, and I 
only. Uncertainties due to the use of only 3 filters have been 
assessed by Monte Carlo simulations (addressed in Section 
[3:311. 



3.4-1 Cluster age and their reliability 

The cluster population of ESO 185 shows a peak at an 
age of 3.5 Myr (log(age)= 6.5; top panel of Figure [9]). The 
strength of this peak is quite robust, since at this age bin 
the light of the cluster is dominated by massive blue stars 
whi ch have short life times and create unique color features 
fsee [Maiz Apellanij (|2009h for more details). We tested our 
models and the uncertainties produced by the fit (see 
the previous section), and none of the combinations of fil- 
ters we used in this analysis would produce a spurious peak 
at 3.5 Myr. In the plot, we point out with vertical lines which 
age bins are affected by uncertainties. A fit performed on a 
combination of the R, I, and H filters produces a spurious 
narrow peak at log(age)= 6.7 (5 Myr; dotted-dashed line), 
and leads to the underestimation of the real ages. However 
the green histogram in the top panel of Figure [51 shows that 
not a single cluster in our sample of the R, I, and H de- 
tections has an estimated age in that age bin. We noticed 
from our simulations that all the used combinations of fil- 
ters produce fake peaks at the age bins log(age)= 7.1 and 
7.3, while emptying the age bins between 30 and 100 Myr 
(from log(age)= 7.5 to 8.0). This effect is due to a loop in 
the color-color space model of the tracks and is commonly 
known to create features of this type. In this age range we 
found 70 objects, so there is a possibility that a fraction of 
them has older ages (up to 100 Myr). This uncertainty is 
taken into account in the further discussions in this work. 
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Figure 9. Final age, mass, and extinction distributions for 264 
clusters in solid black line histograms. The plots include also the 
individual distributions of the 43 clusters detected in R, I, and 
H bands (dotted thick green line) and the 91 clusters for which 
properties were find using FUV, R, and I detections (solid purple 
line). See the main text for details. 



The final age distribution (black continuos line distribution 
in the top panel of Figure [91 shows that the starburst phase 
started less then 100 Myr ago and became very active in the 
last 30 Myr. We also found ~ 10 older (age ^ 1 Gyr) clus- 
ters (globular clusters). Thus, the galaxy is not producing 
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the first generation of stars, but has been rejuvenated by the 
recent merger with a gas-rich dwarf. 



3.4-2 The cluster mass distribution 

The final mass distribution (central plot in Figure [9]) shows 
a peak at a few times 10* Mq. As the simulations show (Fig- 
ure [71 central panels), there is a tendency in our analysis 
to underestimate the mass of some of the clusters. In SED 
fits involving more than 3 filters, the difference is not more 
than 30 %. On the other hand, if only the FUV, R, and 
I filters are used, the masses are possibly underestimated 
up to a factor of 2. The plot in Figure [9] shows that the 
masses estimated from a UV-RI fit have a peak at 10*Mq. 
It is likely that a fraction of those clusters are possibly more 
massive. A counter-tendency is observed in the fit involving 
the R, I, and H bands (Figure [T] central right panel), where 
a third of the low mass cluster population has recovered 
masses overestimated up to a factor of 2, but the effect is 
smaller at higher masses. The most massive clusters in the 
observed mass distribution are determined with a RIH fit 
(green histogram in the central plot of Figure [9} , implying 
the possibility that some of them in reality may be a factor 
of 2 less massive, even if their average mass is more likely to 
be underestimated. 



3.4-3 The cluster extinction distribution 

The distribution of the extinction, E{B — V), peaks close to 
zero and has a slow decline to values higher then l.O (bot- 
tom panel in Figure [9]). We observe that values of extinction 
higher than 0.4 are all constrained with the RIH fit (green 
dotted line distribution). Using simulations (bottom right 
panel in Figure O we found that for this type of fit there 
could be an overestimation of up to a factor of 4 for clusters 
with an input extinction of 0.05 (black hatched histogram). 
However, for higher input values {E{B — V) = 0.2 and 0.3) 
this effect tended to diminish, with overestimations smaller 
than a factor of 2. In Figure [TO] the estimated extinctions 
of the clusters versus the ages are shown. We found that 
many of the clusters with extinctions higher than 0.6 mag 
have ages smaller than 3.5 Myr. At this age range clusters 
are still involved in the partially embedded phase so we ex- 
pect to observe a high scatter in the extinctions. Even if 
some values of E{B — V) are overestimated in our analysis 
(mainly for young clusters detected in R, I and H) we expect 
the uncertainty to be smaller than a factor of 2. 
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Figure 10. Extinction versus age diagram. The blue (grey) aster- 
isks show each single cluster. Superimposed the mean extinction 
estimated at each model age step with the corresponding standard 
deviations. 
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Figure 11. Mass versus age diagram. Clusters which do not show 
any red excess are represented with filled black dots. Clusters af- 
fected by IR excess are showed with red (grey) triangles. The blue 
(grey) diamonds are for clusters affected by I and IR (where de- 
tected) bands excess. The lines (see the inset) show the detection 
boundaries in ages and masses corresponding to the magnitude 
limits reached at ct ^ 0.2 in each filter. As already expected there 
are not clusters sitting below these limits. 



3-4-4 Revealing the properties of the clusters population 

In Figure [TT] we show the masses of the clusters as func- 
tion of the ages. Due to our inclusion of only clusters with 
detections in at least three filter in the SED fit, all clus- 
ters lie above at least 3 detection limit lines. Most of the 
clusters, ^ 60 %, are located at ages below 10 Myr and 
have masses between 10"^ and 10^ Mq. As already observed 
in Haro 11, we did not find low mass (10^ — 10* M©) clus- 
ters at very young ages (below 3 Myr). Two mechanisms, 
blending and extinction (partially embedded phase) can re- 
move objects in this region. Clusters in active star-forming 
environments do not form in isolation but in complexes. At 



the distance of ESO 185 (similarly to Haro 11), it is possible 
that less massive clusters are blended. On the other hand, 
as we observe in Figure 1101 extinction at this age range is 
quite spre ad, indicating that the em bedded phase lasts sev- 
eral Myr (jLarsenI [ioiO : iLadal [ioiol l . The newborn clusters 
need time to clean their surroundings from the natal GMCs 
where they formed. For clusters with masses below 10* Mq, 
even a moderate extinction of E{B — 1/) « 0.5 (Ay ~ 1.6) 
would move the object below the detection limits. In addi- 
tion, a significant fraction of young clusters are affected by 
red excess, in which case the H and/or I filters are removed 
from the SED fit which then has to rely on the shallower (in 
terms of the mass detection limit) U and B bands. Hence, 
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Figure 12. optical (i? — /) versus NIR {R — H) color-color dia- 
gram. The position and the associated photometric error of the 
4 extended and NIR bright cluster candidates are plotted. The 
apparent F184W magnitudes of the clusters # 13, 39, 275, and 
366 are mj = 19.5, 22.5, 20.3 and 21.1, res pectively. The blue solid 
line is the evolutionary track predicted bv lZackrisson et al.l ll200lh 
models at the redshift and mctallicity of the galaxy. Several age 
steps are indicated. The arrows show where the clusters # 39 and 
# 275 will move if a correction for extinction of E{B — V)= 2.12 
and 1.24 mag is applied. The other 2 objects have R — I colors 
impossible to explain with the current models at the redshift of 
ESO 185. 



the apparent dearth of low mass young clusters is probably 

not real. 

Contrary to what found in Haro 11 (|Adamo et al.l 

l2010al ). the most massive clusters are not in the youngest 
bins, but at very old ages (Figure llip . We found signifi- 
cant secondary peaks of clusters with ages around 1 and 14 
Gyr, and masses comparable to the most massive globular 
clusters observed in the Milky Way. Even if the number of 
objects is smaller, the total mass in observed intermediate 
age (0.1-1 Gyr) and old (1-14 Gyr) clusters is comparable 
to that of young clusters, in each case a few times lO^M©. 
However, the mass detection limit is higher for the older 
clusters and the total inferred cluster population mass is 
therefore higher. This suggests that the galaxy has experi- 
enced extreme star formation episodes also in the past. We 
will return to this point and discuss implications for the past 
evolution of ESO 185 in Section 15.31 



3.4-5 Cluster misfits or interlopers? 

We found 4 objects which were undetected in UV-U, only 2 
of them partially detected in B, and with increasing bright- 
ness from R to H band. Since these objects are resolved and 
isolated (in the outskirts of the starburst region) we per- 
formed new photometry (using an aperture of 1.0") and a 
new fit to the SEDs. Two of the sources, ^ 39 and # 275 
(Figure [T^ . are consistent with being deeply embedded star 
clusters. They have fitted reddenings of E{B — V) — 2.12 
and 1.24 (the arrows show how the clusters move in the 
color-color diagram, due to this extinction), ages between 
2 and 3 Myr, and masses of several times 10^ M0, respec- 
tively. The other two sources, # 13 and # 366, have a i? — / 



color, impossible to reconcile with the evolutionary track of 
the galaxy. There is hence a high probability that they are 
background objects (either exceptionally bright ^ ~ 3 Ly- 
man break galaxies or z ~ 0.5 ellipticals) and we excluded 
them from the present analysis. They are not included in 
any of the figures showing our overall results for the cluster 
population in ESO 185. 



4 THE CLUSTER RED EXCESS IN ESO 185: 
DIFFERENCES AND ANALOGIES WITH 
OTHER CASES IN THE LITERATURE 

In this section we will discuss the cause of red excess seen 
in many ESO 185 clusters and compare to other galax- 
ies/environments found in the literature. 

At low metallicities (Z ^ 0.008), the contribution from 
photoionized g as around clusters younger than 10 Myr is 
not negligible i Kriiger et ah 19951: Zac krisson et aL 2001 : 
Bergvall fc Ostlin 2002; Zac krisson et al . 2008; Rcin cs et al.1 



2OIOI : Adamo et al.l l2010blh . If the spectral energy distri- 



butions (SEDs) of the clusters would be fitted by models 
including only stellar continuum, it produces a consider- 
able observed flux excess at NIR bands. The fraction of 
nebular emission contributing to the integrated fluxes in - 
creases and last longer in the NIR f see I Adamo et aill2010bt ). 
The red excess found in the young star clus ters of NGC 
4449 |Reines et al.|[2008al ') and SBS 0335-052E (jReines et al.1 
l2008b ) is pro duced by the nebular emission surrounding the 
clusters (see iReines et al. 2010; Adamo ct al. 2010lJ ). On 
the other hand, the N IR excess found in Haro 1 1 star clusters 
(|Adamo et al.ll2010al ') and ESO 185 sit above models which 
already include nebular contribution so other mechanisms 
need to be invoked. 

In Figure 1111 we differentiate clusters affected by red 
excess using different symbols. With red triangles, we show 
clusters affected by excess in the IR (H band). In blue dia- 
monds we display clusters which possess NIR flux excess, i.e. 
clusters affected by excess at A > 8000 A. Clearly, clusters 
that show a NIR excess are around 5 Myr old or younger, 
and have masses below 10^ M0. Only two clusters affected 
by NIR excess are older (~ 16 — 18 Myr) and more massive 
(> 10^ Mq). Clusters with only IR excess (i.e. with H-band 
but without I-band excess) have a wider range of ages (from 
a few to 35 Myr) and masses (between 10^ and few times 
10^ M0). 

A correlation between the strength of the observed flux 
excess (Am — mmod — n^obs) and the derived ages, masses, 
and/or extinctions can give valuable information on the un- 
derlying mechanisms. In Figure[T3]we show in blue, clusters 
with a NIR excess and in red, those with IR excess. Clusters 
younger than 6 Myr are drawn as filled dots, while the older 
are showed as open squares. The median values of the resid- 
uals in the H and I bands of star clusters not affected by red 
excess are also included for each age, mass, and extinction 
bin (black crosses). 

Studies of very young clusters in the Milky Way show 
that for ages below 6 Myr, the cluster is experiencing a rapid 
and complex evolution from being deeply embedded in the 
natal giant molecular cloud (GMC) to shining in an envi- 
ron ment that has been cleared from dust and molecular gas 
(see iLada fc Ladall2003l for a review) . During this period the 
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Figure 13. Strength of the flux excess in If and f bands as func- 
tion of the cluster properties derived from fitting the blue side of 
the SEDs. The vertical dotted lines in the left side panels sepa- 
rate clusters younger than 6 Myr (filled dots) from older clusters 
(open squares). In red we show clusters affected by IR excess (H 
band), and in blue the ones with NIR excess (I and H bands). 
The crosses show, for each bin, the median value of the residuals 
for normal clusters in the sample. 



short-lived massive stars emit strong UV-radiation which 
ionise the ambient gas and heat the dust. The presence of 
dust in the cluster medium and strong UV-fields can cre- 
ate a rise of the continuum at two difTerent wavelengths. 
Hot dust usually re-emits at wavelengths much longer than 
1 ^m, of course sensitively depending on the dust tempera- 
ture. At shorter wavelengths, between 0.6 — 0.9 fj,m a pho- 
toluminescence process involving dust grains and hard UV- 
fields can create the so-called extended red emission (ERE, 
IWitt fc Viih ,2004.). ERE as been observed as a soft " rise" in 
the continuum at 0.6—0.9 / xm of Hll regions and gaseous neb- 
ula ( among many others Gordon et al. I I2OOOI : [ Perrin et al.l 



i (a , 

ll995l : |Perrin fc SivanHiggi T 

Moreover, in a newborn cluster a large fraction of stars 
may still be in the contracting phase as pre-main sequence 
(PMS) objects. iMever fc GreissJ (|2005l ) estimated that PMS 
stars contributed between 1 and 17 % of the integrated H 
band flux of an unresolved cluster, depending of the stel- 
lar popul ation age (from 1 t o 10 Myr) and the used IMF. 
Recently. ICreissl et all (|2010l ) found direct evidence of PMS 
stars in NIR spectroscopic data of a young knot (89/90) in 
the Antennae system. Among the pre-main sequence stars, 
the young stellar objects (YSOs) show a strong NIR excess 
caused by their luminous accretion disks. Interestingly, a 
large number of YSOs are detected at the edges o f young star 
form ing regions, due to triggered star formation (|Elmegreenl 
1 19981 ). These transition phases are typically o bserved in 
mass i ve star-forming regi ons of the Milky Way (|Bik et al.l 
I2OIOI : iMartins et al. I2OI0I ) and very yo ung star clusters in 
the Small Mage Uanic Cloud (NGC 602 , ICarlson et al.|[2007l : 
NGC 346/N66. [Gouliermis et allboiol ). 

Resolved nearby clusters give further insights into the 
main sources of the NIR excess: ERE, at shorter wavelengths 
(~ 0.8^m), hot dust and a large fraction of YSOs and PMS 
stars in the IR (> 1.5/im) . The young knots of the nearby 
starburst spiral NGC 253 (|Fernandez-Ontiveros et al.ll2009l ) 
show a soft bump at 1 — 2 nm which could be reproduced by 



an substantial fraction of YSOs. Studies of embedded clus- 
ters in nearby dwarf starburst galaxies have also revealed 
IR excess w hich has been int erpreted as d ue to hot dust 
(NGC 5253, ICresci et al.]|2005l: Henize 2-10. ICabanac et all 
I2OO5I ). Clusters in Haroll jAdamo et al.]|2010al) . and pos- 
sibly also in SBS 0335-052 l| Adamo et ai1l2010bl ). show an 
I-band excess which may be explained by the ERE feature. 

In Figure 1131 clusters represented by dots (ages below 
6 Myr) are the ones which potentially could be affected by 
these processes. The strength of the excess in the H band 
appears to be higher for younger, lower mass and less extinct 
clusters. However the scatter is large (top panels) and there 
is a general decrease of the number of clusters with increas- 
ing age, mass and extinction (see Figure [9]). The conclusion 
is that H-band excess seems to affect an approximately ran- 
dom (with respect to mass and extinction) 10 % fraction of 
young (age ^ 30Myr) clusters . If the excess is caused by hot 
dust in the cluster interstellar medium we would expect to 
find higher foreground extinction, resulting from a clumpy 
distribution in front of the cluster, shielding partially the 
light produced by the young stars. If there is no/little dust 
along the line of sight to the clusters (very low extinction), 
it is likely that other sources are dominating and the best 
candidates are YSOs. The aperture used to do photometry 
(radius 35 pc) is much larger that the expected radius of the 
clusters (5 — 10 pc). Therefore, inside the region where we 
do photometry, we are including not only the cluster, but 
potentially also star formation at its rim, or po ssibly other 
newb orn clusters, since they form in complexes jElmegreenl 
I2OIOI ). The edges of very young clusters are places for trig- 
gered star formation (large fraction of YSOs) and these re- 
gions may be cleared from their dust at a time when the cen- 
tral massive cluster has already been naked for some time. 
So, there is a fair chance that the excess we are observing 
in the H band is possibly a combination of YSOs and hot 
dust. The bottom panels show the I band excess as function 
of ages, masses and extinctions. No particular trend is ob- 
served, except that the excess at this waveband is only seen 
for very young star clusters, showing that the ERE effect is 
probably the best explanation for this behavior. 

At older ages (> 7 Myr) red super giants (RSGs) af- 
fect considerably the integrated flux es at wavelengths red- 
der than 10000 A (|Lancon et al.ll2008l ). They may be relevant 
for the red excess for two reasons: Firstly, at low metallici- 
ties, the models have a tendency to underpredict the number 
of RSGs observed and there is evi dence that their number 
incr eases at lower metallicity (see lEggenberger et al] I2OO2I 
and iMaeder fc Mevnetl I2OOII for details). Since the RSGs 
dominate the hght in the NIR, if the models underestimate 
their number the observed integrated fluxes would shown 
an excess at such wavelengths. Secondly, in less massive 
clusters (< 10^ M^) stochas tic effects in the IMF sampling 
(ICervino fc Luridianall2004 ) will severely affect the colours 
of clusters, in particular in the IR. Only clusters with mass 
on the order of 10^ Af© can fully sample the IMF[f|- Such sta- 
tistical fluctuations are likely to occur and affect the NIR 
SED when the expected number of RSGs is on the order of 
unity or less. In Figure [T31 clusters plotted as squares (7-35 



^ e.g. for a Salpeter IMF from 0.1 to 100 Mq there will be one 
star more massive than 80 A/q per 10^ A/0 of cluster mass 
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Myr) are the ones where RSGs could have an impact on the 
flux excess. The masses of some of these clusters are low 
enough that IMF sampling effects may be important, and 
the varying number of RSGs from cluster to cluster could 
explain the spread in Ampieow ■ 



5 THE STARBURST PROPERTIES OF ESO 
185 AS REVEALED BY THE STAR 
CLUSTERS 

5.1 The cluster luminosity function 

The cluster luminosity function (CLF) is a useful tool to 
study the formation and evolution of clusters in galaxies. It 
has been shown for many nearby galaxies that the cluster 
mass function (CMF) can be approximated by a power-law 
with index appro ximately a ~ —2 . (M 51, feik et aLll2003l : 
the Antennae, see IWhitmore et al.ll201(]| for a new revisited 
CLF and C MF of the system; LM C and SMC. lHunter et al.1 
I2OO3I : M83, IChandar et"alll2010l '). Assuming that the CLF 
is built up of fully sampled CMFs formed at different ages, 
then the resulting CLF is supposed to have the same power- 
law shape. In reality the observed CLFs have a large spread 
of power law indexes ranging from —2.4 q ^ —1.7 (among 
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have a double power law CLF, steeper than —2.0 at the 
bright end, and approximately —2.0 at the faint end. The 
bend occurs between roughly —8.0 and —10.0 mag (Gieles 
et al. 2006a,b). The steepening of the CLF (both for single 
power laws and for the bright-end of the double ones ) could 
be ca used by a truncation at high mass in the CMF (|Larsenl 
120021 '). or in other words, a limit on the maximum mass 
with which a cluster can be formed. ICieles et all (|2006bl ') 
showed that galaxies whith a double power law CLF im- 
pose physical limits on the formed maximum cluster mass. 
A steepening of the CLF as function of the wavelengths is 
expected if there is a truncation in the CMF, since evolved 
stellar popula tions are brighter at red wavelengths. Recently, 
iLarsenl (|2009l ) suggested that instead of a sharp trunca- 
tion in the CMF, a Schechter function with a characteristic 
mass, Mc « 2 x 10^ M0, and an index of —2.0 at the low- 
mass end is a better representation of the CMF in quiescent 
spiral galaxies. Shear and/or regular patterns (like spiral 
arms) effectively limit the mass of GMCs from which clus- 
ters form. In presence of merging events which enhance the 
compression and the star formation efficiency, the Mc has 
likely higher values. A Schechter-like CMF could explain the 
steepening of the CLF at the bright en ds, the possibility for 
galaxies to form ver y massive clusters (|Bastian et al.|[2006l : 
lAdamo et al.|[2010al ). and why the most massive young clus- 
ters scal e with the glo bal star formation efficiency of the 
systems (|Bastianll2008l ). 

ESO 185 has produced many clusters in the most recent, 
and still active, burst episode. Using the two deep WFPC2 
exposures we detected more than 400 clusters (see Section 
I2.2.1|) . Of these, we were able to derive ages, masses, and 
extinctions of roughly 50 %. In Figure [14] we show the CLFs 
recovered for the F814W and the F606W filters for the en- 
tire population of clusters detected in both bands. The com- 
pleteness tests (Section l2.2.2p show that in crowded regions. 
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Figure 14. Cluster luminosity functions in the two deepest 
frames available for ESO 185 (R and I filters). The dashed his- 
tograms are the observed distributions, the solid ones have been 
corrected for completeness. The power law fit performed on bins 
with 100 % completeness is shown in black. In green (grey) we 
show the fit to the luminosity bins brighter than -11.0 mag. 



already at 26.0 mag a considerable fraction of objects is not 
recovered. Using a conservative limit we fitted the CLF us- 
ing magnitude bins brighter than 25.0 mag in both filters 
(the bins have been corrected for completeness, black con- 
tinuum line histogram) . This first fit is showed with a thick 
black line in the two panels of Figure [141 The power-law in- 
dexes found are indicated and consistent with a « —1.8 in 
both filters. It is clear that the number of objects decrease 
rapidly in bins fainter than 25.0 mag. We checked whether 
including the total number of clusters detected in both fil- 
ters (1265 objects) would change the slope at the faint ends 
(magnitude bins > 25.0 mag). We recovered a quite fiat 
slope of ~ —1.2 at such faint magnitude bins, while the 
slope at brighter magnitudes was unchanged. Finally, we 
tried a fit to the very bright magnitude bins (brighter than 

— 11.0 in absolute magnitudes). Bins brighter than —11.0 
mag are usually scarcely populated. The brightest clusters 
in th e Antennae system are ~ —11.5 mag ([Whitmore et al] 
[2010t ): nearby spiral galaxi es, like M83 and M5 1 , barely fill 
the — 11.0 magnitud e bin ([Chandar et al.l [2OIOI : [Haas et al.l 
I2008h . l6stlinl(|2000l ) showed that blue compact galaxies are, 
instead, very efficient in producing clusters brighter than 

— 11.0 mag. HST data have revealed that small starburst 
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galaxies with total stellar masses 10^ — 1 0^° Mp) have a large 
fraction of very ma ssive young clusters (|6stlin et alj 12003 : 
lAd amo et al ] |2010al : this work). The fit to these bright bins, 
in ESO 185, produced a slope a factor of 0.5 steeper (green 
dotted line). Due to the distance of the galaxy (~ 80 Mpc), 
even the high resolution power of the PC camera, it's not 
enough to spatially resolve all the single clusters. We expect, 
in fact, that not all the clusters form in isolation. Blend- 
ing and crowding are affecting our data analysis. Usually, 
these mechanisms tend to flatten the slope at the faint-ends 
with a few percent, while the slope at the bright bins is 
almost unchanged. Fully understanding the observed lumi- 
nosity distribution of the clusters in ESO 185 (flat at low 
luminosities and approximately —2.0 at bright ends) is chal- 
lenging. Monte Carlo simulations are needed to investigate 
in more detail the possible causes (Adamo et al. in prep.). 
We discuss, here briefly, possible scenarios. 

Blending, a cha nge of slope in the CMF at low masses 
l|Gieles et al.ll2006bh , and a mass dependent cluster disrup- 
tion timescale (low mass clus ters dissolve faster than high 
mass ones; see iLamerj [ioogi for a short review) could all 
produce a flattening of the CLE. The mass-age diagram 
in Figure 1111 shows that there are empty regions at low 
mass (^ IO^Mq) at very young ages (1 — 3 Myr) and be- 
tween 20 Myr and 100 Myr. We already discussed that 
very young low mass objects could be blended or fall below 
the detection threshold if affected by moderate extinction 
{E{B - V) > 0.5). At older ages, between 30-100 Myr, we 
expect that we might be seeing to few clusters since many 
of them have been assigned younger ages (~ 20 Myr) by the 
SED flt (see Section ??). However, the CLF is not affected 
by systematics in the performed cluster age analysis and rep- 
resents an independent test of the cluster formation. If we 
introduced important biases in the age estimates this would 
be revealed by a normal recovered CLF with index —2.0. 
However, we observe an important flattening of the CLF 
which suggests that the missing fraction of low mass clus- 
ters are either heavily blended or preferentially destroyed. 
None of the two scenarios can be excluded. F inally, we al- 
ready noti ced in Haro 11 jAdamo et al.l l 2010al ) and in ESO 
338-IG04 (jOsthn et ai]l2003l l that the CLF in blue compact 
galaxies tend to be flatter than the typical values (found in 
spirals and interacting systems) and produce high fractions 
of massive clusters. It's not possible to draw any conclusion 
from the available data, but it appears that the galactic en- 
vironment in these dwarf starburst galaxies has favored the 
formation of massive clusters and/or the disruption of low 
mass clusters, playing a crucial role in shaping their CLF. 

5.2 CFR versus SFR in ESO 185 



Using the total Ha luminosi tv (lOstlin et al.lll999l) and the 
radius of ESO 185 (3.8 Kpc. lOstlin et aLlbOOlh Twe applied 
the Kennicutt relation (jKennicuttl 1 1998f l to obtain an esti- 
mate of the present global star formation rate in the galaxy: 
SFR« 6.4 ± l.OMgyr"^; and of the SFR surface density, 
EsFR ~ O.14M0yr~'^pc~^. Assuming that a whole cluster 
population has formed every 10 Myr and that we are com- 
plete in detecting all clusters more massive than 10* Mq, we 
derived the observed cluster formation rate, CFR. We as- 
sumed that the clusters have been formed with a power law 
CMF of index —2.0 and estimated the total CFR by extrapo- 



lating to clus ters less massive than 10*' Mq (a similar exercise 
was done in lAdamo et all l2010al following ICoddard et al.l 
l2010h . 

Since we know the present SFR and CFR in th e galaxy, 
we ca n estimate the cluster formation efficiency, F (jBastianl 
|2008| ). F is defined as the ratio between the CFR and SFR 
and is an estimate of the fraction of star formation hap- 
pening in compact clusters. For ESO 185 we derived a 
F = 0.26±0.05, or in other words, 26±5 % of the star forma- 
tion is happening in clusters. Although this r atio is smaller 
than t he one found for Haro 11 (see Table 4 in I Adamo et al.l 
l2010al ) , is still high compared to many ot her starburst sys- 
tems. Interestingly, ICoddard et al.l (|2010l ) showed that the 
star formation density of the host systems and the cor- 
responding values of F are positively correlated, i.e. that 
higher star formation rates in the host galaxy enable the 
production of more numerous and more massive clusters. 
As already done for Haro 11, we checked whether ESO 185 
followed the same relation. In the right panel of Figure [151 
we see that the position of ESO 185 is only 2a above the 
exp ected value p r edicte d by the relation (18 %). 

lOsthn et al.] (|200ll l produced an approximative esti- 
mate of the burst mass (i.e. the mass in young stars) in 
ESO 185: 2.ltofs, x 10* Mq. If we consider the extrapolated 
total mass in clusters younger than 40 Myr (~ 3.4 x lO'' Mq) 
we estimate that roughly 16 % of the burst mass has been 
produced in bound clusters. 

In Figure [151 we also show the cluster formation history 
in the last 40 Myr of ESO 185. The uncertainties in the esti- 
mation of the cluster ages between 20 and 10 Myr make the 
rate of cluster formation at those ranges not fully reliable. 
It is clear, however, that the galaxy is experiencing now a 
peak of cluster formation (l.TMoyr^^). 



5.3 ESO 185 progenitors and the starburst 
evolution 

The analysis of ESO 185 morphology and gas dynamics re- 
veals a galaxy that has experienced a merger event. The 
recent cluster production in the galaxy was probably caused 
by the infall of gas in the merger event. The discovery of a 
shell structure in the deep WFPC2 images and the numer- 
ous clusters in the galaxy could be used to trace the possible 
progenitor s and evolution of t he system. From lOstlin et ahl 
pool.) and lOstlin et al. I (|l999t ) we know that a second com- 
ponent in the galaxy is counter-rotating and has a mass 
which is ~ 1/10 of the main component. Since the total 
stellar mass of the system is 7 x 10^ Af© we expect the the 
merger process has involv ed dwarf systems. Even th ough the 
simulations performed bv lKoiima fc Noguchil (|l997l ) did not 
consider merger between dwarfs, we make a qualitative com- 
parison with their results. From their 4 different sets of sim- 
ulations, the ones produced by radial encounters are most 
similar to the morphology of ESO 185; the merging disk 
produces a shell formed of both stars (from the early type 
galaxy) and gas (from the sinking satellite), while a tidal 
stream is left behind and a bar-like structure is formed in 
the center of the potential. The simulation with a smoothed 
potential well shows short-lived shell structures. They no- 
ticed that when the shells form the star formation in the 
galaxy is ceased due to the dissipative nature of the gas 
particles. In the real case of ESO 185, we derived the star 
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Figure 15. Left panel: Cluster formation rates during the last 40 Myr of starburst activity. The filled dots connected by the green 
dotted line show the observed CFR derived from clusters more massive than 10'' M0. The thick squares show the derived CFR, if the 
total mass in clusters less massive than 10* is extrapolated using a CMF with i ndex —2.0 down to 10 '^ Mp;. Right panel: The cluster 
formation efficiency (F) vs the star formation rate surface density, using data from lCoddard et al. with the positions of Haro 11 

l|Adamo et al.ll2010a^ and ESO 185 added. Sec text for details. 



formation history in the galaxy, following the formation of 
clusters. In Figure [161 we show how clusters of different age 
ranges are located in the galaxy. The most evolved clus- 
ters, the globular clusters, older than 0.1 Gyr, appear to be 
biased towards the south part of the galaxy. However this 
could be a projection effect, since we don't know the incli- 
nation of the galaxy. Two of these clusters are in the shell 
and have ages of 6 and 14 Gyr. These globular clusters were 
most likely part of the early type progenitor and in the oscil- 
lations around the potential that formed the shell dragged 
together with the stars. The panel in the bottom left side 
show clusters in the age range 20-100 Myr. A fraction of 
objects from the top right panels should in reality be in this 
older range. We observed that the clusters at 20-100 Myr 
and at 10-20 Myr age ranges are more central distributed. 
Two of the clusters in the shell with detection in more than 
2 filters, have ages of 15-20 Myr. Clusters between the shell 
and the starburst region are younger than 10 Myr and visible 
in the top left panel of Figure [161 This panel shows also that 
in the last 10 Myr the starburst has rapidly expanded out- 
ward from the central region. The final encounter between 
the galaxies happened between 100 and 20 Myr ago. The 
high rate at which the clusters have been formed in the last 
20 Myr proves this scenario. In fact, the merging event has 
probably enhanced the pressure in the intergalactic medium 
and the compression of the gas, making the star formation 
more efficient. The shell is at a projected distance of only 5 
Kpc from the center of the galaxy. Such proximity suggests 
that it has recently formed. Concentric with the shell and 
at approximately one third of the distance from the centre 
another bent narrow structure is visible that could possibly 
be an inner shell, or possibly a small spiral arm. If we as- 
sume that the outer shell (marked in Figure 1) was the first 
one and note that young clusters are found only out to the 
shell, b ut not beyond it, the shel l must be no older than 
10 Myr. IKo] ima fc Noguchil (|l997l ) noticed that the forma- 
tion of the shell happens when the system is entering in a 
po st-starburst phas e. Simi lar conclusions have been reached 
bv iMavbhate et all (|2007l l. who observed that the ages of 
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Figure 16. Position of the clusters in the galaxy as function 
of their ages. Clusters affected by red excess are showed in red 
triangles (excess in H band), and in blue diamonds (excess in I 
and H bands). The center of the galaxy and the position of the 
shell are shown by a green cross and arc. 



the globular clusters in the shell correspond to the youngest 
stellar population formed in the galaxy before it went into 
a more quenching phase. It is plausible that after this burst 
episode, ESO 185 will experience a more quiet post starburst 
phase. However, the galaxy has a mass of neutral hydrogen, 
Mh7 ~ 3 X 10^ M0 (BergvaU et al., in prep.). At the current 
SFR (?a 6.4 ± l.OMoyr"^), the present burst could be ac- 
tive for another 500 Myr, although the fraction of gas close 
to the centre and currently available for star formation is 
unknown. 

While the cluster formation in ESO 185 is now very ac- 
tive the most massive clusters that we find in this galaxy 
are found among the > 100 Myr old clusters. The total 
mass in clusters younger than 100 Myr and more massive 
than a mass completeness limit of 10* is 2 x IO'^Mq. As- 
suming a power- law mass function with index a = — 2 for 
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masses below we estimate a total mass in young clusters of 
3.4 X 10^ Mq. We note that no mass has been ascribed to 
the red excess luminosity. 

In the age range 100 Myr to 1 Gyr we find a total cluster 
mass of 1.9 x IO^Mq when counting clusters more massive 
than IQ^Mq. If we assume that clusters in this age range 
has a Gaussian CMF similar to that of Milky Way globular 
clusters (Mp^ak = 2 x 10^ Mq, a = 0.6 log(M) [Mq]) the 
expected mass in clusters below 10^ is negligible. In this age 
range, and as previously noted, we see what appears to be a 
dearth of clusters with mass in the range a few times IC* to 
10^, but this has little impact on the total estimated mass of 
this sub-population (unless there are in reality many more 
low mass clusters than seen and accounted for). 

For old (age > IGyr) clusters with properties sim- 
ilar to GCs we see a total mass of 3.1 x lO^Afc when 
including clusters with masses down to 5 x 10^ A/©. As- 
suming the same Gaussian CMF as in the previous exam- 
ple, the total estimated mass for the old GC population is 
1.1 x 10*^1/0. The stellar mass of the g alaxy underly i ng the 
starburst component was estimated bv lOstlin et al.l (|200lh 
to Mgai « 6 X 10^ Mq which would imply that the GCs in 
ESO 185 make up more than a percent of the total stellar 
mass. This value is high in comparison with most galaxies, 
and a similar result was found for ESO 338-04 (|6stlin et all 
I2OO3I ) and may indicate a greater survival rate for GCs in 
low mass systems. If we instead of the cluster population 
mass view the number of old clusters, we get a similar pic- 
ture. The total number of old clusters is 11. Correcting for 
completeness we expect 30 (again assuming a gaussian LF 
identical to that for the M W). From this we find a mass nor- 
malised specific frequency (| Ashman fc Zepj|l998h of T f« 5, 
which is typical for early type systems, and consistent with 
the picture that ESO 185 as we see it today is the product 
of a merger between an early type low mass galaxy and a 
gas-rich dwarf. 



6 CONCLUSIONS 

This work, togethe r with the analysis of ESO 338 -IG04 
lOstlin et all |2003| ). Haro 11 (|Adamo et al.1 bOlOal ). and 
MRK 930 (Adamo et al. in prep.) provides a statistical study 
of the properties of star clusters in blue compact galaxies 
and of their formation and evolution in general. 

In this paper we have used images obtained with the 
HST in several passbands ranging from the UV Xcen = 
0.15/im to the near IR (1.6/^m) of the luminous blue compact 
galaxy ESO 185. The new HST images confirms the presence 
of a tidal tail/plume and in addition finds a shell, corrobo- 
rating that this starburst was triggered by a (dwarf) galaxy 
merger. We find several hundreds of luminous star clusters 
which we analyse by means of SED fitting using state of the 
art models including a nebular gas emission. The clusters 
span a wide range of ages, from newborn < IMyr to old 
(lOGyr) globular clusters, and the masses range between a 
few times 10'^ to IO^Mq. A recent peak in the cluster forma- 
tion is visible at an age of 3.5 Myr, and at least 20 % of the 
present star formation is happening in compact (unresolved) 
star clusters. 

The dust extinction for young clusters present a large 
scatter, but rapidly decrease for ages > 5Myr. This indicates 



that the embedded phase does not last for a fixed la pse of 
time, but is a tran sient phase in the cluster evolution (|La di 
l20ld : lLarsenll2010l '). 

We found that a fraction of clusters younger than 6 
Myr displays a flux excess at A > 8000 A, so-called NIR ex- 
cess. We discuss possible scenarios which could explain why 
the observed fluxes sit above the prediction of evolution- 
ary models. Hot dust and a considerable contributions from 
young stellar objects (YSOs) are among the advocated ex- 
planations for similar excess cases found in the literature. 
Resolved stellar populations of young star clusters in the 
SMC and in the Milky Way, show that even if the cen- 
tral part of the system is gas and dust free, star forma- 
tion, t riggered by posit ive feedback, is still active at the 
edges (jElmegreenl Il998l ). Considering that the very young 
star clusters in ESO 185 are quite massive and that the pho- 
tometric aperture used is spatially bigger than the real size 
of the compact clusters, we are likely including contribu- 
tions from the edges of the star forming regions where hot 
dust and YSOs reside in our photometry. The clusters in 
ESO 185, for which we observe a NIR excess, may therefore 
be objects that are still in a partially embedded phase. This 
scenario could also explain why we detect those clusters at 
optical wavelengths. 

Another fraction of clusters with ages older than 6 Myr 
showed a flux excess only for wavelengths ~ 1.6 ^m, referred 
to as IR excess. In this case, we conjecture that Red Super 
Giants (RSGs) may contribute to the excess. 

The cluster luminosity function (CLF) becomes steeper 
at brighter luminosity and is very flat at the faint ends. A 
single power law fit gives an index of —1.8, which however is 
a rather poor representation of the CLF. We discuss possi- 
ble mechanisms which could cause a flattening of the CLF. 
Blending and projection effects are addressed, together with 
a possible missing fraction of low mass clusters at very young 
ages (below 3 Myr) and around 30 Myr and 100 Myr. 

The tidal stream and shell characterize the morphol- 
ogy of ESO 185, and help us to trace back the progenitors 
of the pres ent syste m. U sing hydrodynamical simulations, 
iKoiima fc Noguchi (| 19971 ') showed that a minor merger be- 
tween an elliptical and a gas rich satellite, create shells of 
stars and gas, and tidal tails. Based on the age of the star 
clusters, the final encounter likely happened between 100-20 
Myr ago. The clusters in and close to the shell are not older 
than 20 Myr. The shell is approximatively 5 Kpc away from 
the cen ter of the galaxy, proving that it has been recently 
formed. iKoiima fc Noguchil (|l99^ observed that the forma- 
tion of the shell and the ceasing of the star formation were 
coeval. Similarly, iMavbhate et al.l (|2007l ) found the younger 
globular clusters in the post-starburst galaxy AM 0139-655 
in the shell. ESO 185 will likely soon transit into a post burst 
phase. There is, however, a considerable fraction of HI that, 
at the current SFR, could support the starburst for the com- 
ing 500 Myr. 
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Figure 1. Continue from Figure [5] 
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Figure 2. Continue from Figure [5] 



Table 1. Continue from Table |2] 



IDs 


Xr 


Q 


Myr 


10^5 Mq 


E{B-V) 


UV-UBRIHfit 


56 


4.6 


0.008 


19.5 


1.4 


0.04 


90 


1.4 


0.392 


1.5 


0.2 


0.08 


99 


7.6 


0.002 


8.5 


0.4 


0.0 


118 


5.1 


0.004 


19.5 


2.0 


0.06 


141 


0.6 


0.727 


9.5 


1.5 


0.05 


146 


2.5 


0.115 


1.5 


0.5 


0.08 


147 


0.6 


0.736 


3.5 


0.4 


0.02 


216 


2.9 


0.066 


19.5 


2.0 


0.06 


246 


3.8 


0.055 


9.5 


0.9 


0.01 


338 


2.9 


0.1205 


4.5 


0.3 


0.07 


UV-UBRI&t 


56 


5.8 


0.009 


19.5 


1.4 


0.04 


90 


1.3 


0.460 


1.5 


0.1 


0.08 


99 


7.6 


0.002 


8.5 


0.4 


0.0 


118 


6.6 


0.004 


19.5 


2.0 


0.06 


141 


0.6 


0.727 


9.5 


1.5 


0.05 


146 


3.6 


0.063 


1.5 


0.6 


0.09 


147 


0.6 


0.736 


3.5 


0.4 


0.02 


216 


2.6 


0.160 


19.5 


1.9 


0.05 


246 


3.8 


0.055 


9.5 


0.9 


0.01 


338 


2.9 


0.1205 


4.5 


0.3 


0.07 


UV-UBRIHfit 


56 


3.7 


0.159 


3.5 


0.2 


0.07 


90 


1.8 


0.340 


1.0 


0.2 


0.06 


99 


1.8 


0.410 


3.5 


0.1 


0.02 


118 


3.7 


0.160 


3.5 


0.2 


0.09 


141 


1.3 


0.524 


9.5 


1.5 


0.05 


146 


3.8 


0.152 


1.0 


0.5 


0.06 


147 


0.8 


0.654 


4.5 


0.4 


0.0 


216 


4.0 


0.132 


18.5 


1.6 


0.04 


246 


1.0 


0.601 


3.5 


0.3 


0.06 


338 


1.7 


0.422 


2.5 


0.2 


0.06 
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Table 2. Photometric properties of the clusters showed in Figure [5l[T] and [2] Magnitudes are in Vega system throughput. The kind of 
excess found in the cluster is indicated in column 2. The errors gave in the table includes photometric errors and uncertainties due to 
aperture correction. 



IDs Obs Excess F140LP F336W F439W F606W F814W F160W 



56 


I, H 


21, 


,10 


± 


0, 


,09 


22 


,00 




0, 


10 


23.52 


± 


0.14 


23.29 ± 


0.11 


22 


,96 


± 


0, 


,06 


21 


.40 


± 


0, 


,38 


90 


H 


21, 


,51 


± 


0, 


,09 


22, 


,88 


± 


0, 


19 


23.87 


± 


0.16 


23.07 ± 


0.11 


23, 


,78 


± 


0, 


,13 


22, 


.55 


± 


0, 


,42 


99 


I 


20, 


,70 


± 


0, 


,10 


21, 


,92 


± 


0, 


11 


23.52 


± 


0.15 


23.19 ± 


0.12 


22, 


,75 


± 


0, 


,07 












118 


I, H 


20, 


,93 


± 


0, 


,11 


21, 


,76 


± 


0, 


,09 


23.31 


± 


0.16 


22.99 ± 


0.12 


22, 


,56 


± 


0, 


,09 


20, 


.93 


± 


0, 


,41 


141 


None 


20, 


,18 


± 


0, 


,08 


21, 


,48 


± 


0, 


10 


22.41 


± 


0.08 


22.01 ± 


0.10 


21, 


,63 


± 


0, 


,05 












146 


None 


20, 


,17 


± 


0, 


,06 


21, 


,49 


± 


0, 


,08 


22.59 


± 


0.08 


21.73 ± 


0.10 


22, 


,44 


± 


0, 


,05 


22 


.04 


± 


0, 


,39 


147 


None 


19, 


,75 


± 


0, 


,06 


21, 


,08 




0, 


,08 


22.32 


± 


0.09 


22.28 ± 


0.11 


22, 


,55 


± 


0, 


,07 












212 


I, H 


19, 


,66 


± 


0, 


,06 


20, 


,76 


± 


0, 


,07 


22.21 


± 


0.08 


21.93± 


0.10 


21, 


,61 


± 


0, 


,05 


20, 


.37 


± 


0, 


,39 


215 


H 


21, 


,38 


± 


0, 


,18 


21, 


,87 


± 


0, 


,11 


23.20 


± 


0.10 


22.65 ± 


0.10 


22, 


,80 


± 


0, 


,07 


21, 


.04 


± 


0, 


,38 


216 


H 


21, 


,01 


± 


0, 


,20 


21, 


,84 


± 


0, 


12 


23.01 


± 


0.11 


23.18 ± 


0.16 


22, 


,61 


± 


0, 


,14 


20 


.44 


± 


0, 


,41 


246 


I 


20, 


,35 


± 


0, 


,11 


21, 


,49 


± 


0, 


11 


22.90 


± 


0.15 


22.48 ± 


0.12 


22, 


,10 


± 


0, 


,10 












286 


None 


18, 


,98 




0, 


,05 


20, 


,62 


± 


0. 


,06 


21.63 


± 


0.06 


21.26 ± 


0.09 


20, 


,90 


± 


0, 


,04 


19, 


.92 


± 


0, 


,39 


338 


I 


21, 


,02 


± 


0, 


,14 


22, 


,12 


± 


0, 


,12 


23.60 


± 


0.19 


23.04 ± 


0.11 


23, 


,07 


± 


0, 


,10 












344 


I 


20, 


,41 


± 


0, 


,10 


21, 


,56 


± 


0, 


11 


23.05± 0.16 


22.51± 


0.14 


22, 


,35 


± 


0, 


,10 
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